Branching Brownian Motion 
with catalytic branching at the origin 

Sergey Bocharov and Simon C. Harris 
University of Bath 

February 19, 2013 

Abstract 

We consider a branching Brownian motion in which binary fission takes place only when 
particles are at the origin at a rate /3 > on the local time scale. We obtain results regarding 
the asymptotic behaviour of the number of particles above Xt at time t, for A > 0. As a 
corollary, we establish the almost sure asymptotic speed of the rightmost particle. We also 
prove a Strong Law of Large Numbers for this catalytic branching Brownian motion. 

1 Introduction 

1.1 Model 

In this article we study a branching Brownian motion in which binary fission takes place at the 
origin at rate f3 > on the local time scale. That is, if (X t : t < t) is the path and (L s : s < r) 
is the local time at the origin of the initial Brownian particle particle up until the first fission 
time r, then the first birth occurs at the origin as soon as an independent exponential amount 

of local time has been accumulated with L T = Exp(/3) and X T = 0. Once born, particles 
move off independently from their birth position (at the origin), replicating the behaviour of the 
parent, and so on. Heuristically, we have an inhomogeneous branching Brownian motion with 
instantaneous branching rate /3(x) :— f3 Sq{x) , since we can informally think of Brownian local 
time at the origin as L t = L So(X s ) ds, where Sq is the unit Dirac-mass at 0. 

Although BBM models have been very widely studied, the degenerate nature of such catalytic 
branching at the origin means that the above BBM model needs some special treatment. Related 
models with catalytic branching have been extensively studied in the context of superprocesses; 
for example, see Dawson & Fleischmann [3], Fleischmann & Le Gall |9 or Englander & Turaev [7J. 
In the discrete setting, catalytic branching random walk models have recently been considered 
by, for example, Carmona & Hu [3J and Doring & Roberts [5]. 

1.2 Main results 

In this section, after first setting up some notation, we will state our main results for BBM with 
catalytic branching at the origin (presenting them in the order that we will prove them). 

We denote the set of particles present in the system at time t by N t , labelling particles 
according to the usual Ulam-Harris convention. If u G N t then the position of particle u at time 
t is Xf and its historical path up to time t is (A")o< s <t. Also, we denote the local time process 
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of a particle u £ N t by (L")o< s <t- The law of the branching process started with a single initial 
particle at x is denoted by P x with the corresponding expectation E x . 
Firstly, we shall calculate the expected population growth. 

Lemma 1.1 (Expected total population growth). For t > 0, 



E^Ntlj = 2$(^Vt)e^* - 2e^ 



• 1 as t 



Lemma 1.2 (Expected population growth rates). For A > 0, let Nf* :— {u £ N t : A" > Xt} be 

the set of particles that have an average velocity greater than X at time t. Then, as t — > oo, 

i^z^^A ^ A f \P 2 -HX tfX</3 



-logE(\N t xt \) 



I \2 



X 2 ifX>(3 



Note, the expected growth rate of particles with velocities greater than A > 0, A\, is positive 
or negative according to whether A is less than or greater than (3/2, respectively. That is, the 
expectation speed of the rightmost particle is /3/2. (Also note, by symmetry, similar results hold 
throughout for particles with negative velocities.) 

Next, we consider the almost sure asymptotic behaviour of the population. 

Theorem 1.3 (Almost sure total population growth rate). 

lira = —p F-a.s. 

t— >-oo t 2 

Theorem 1.4 (Almost sure population growth rates). Let X > 0, then: 
(i) ifX>§ then lim^^ \N^\ — OP- a.s. 



log |jV. A *l 

(ii) ifX < § then lim^ * = A A = \p 2 - /3X P-a.s. 

From Theorem 11.41 we immediately recover the speed of the rightmost particle, 



Rt := sup A t " , t > 0. 
Corollary 1.5 (Rightmost particle speed). 

v Rt P p 

lim — = — F-a.s. 

t— s-oo t 2 

We can also say something about the rare events of |A t A '| being positive when we typically 
do not find particles with speeds A > 4. 



Lemma 1.6 (Unusually fast particles). For A > 



lim lo g P(l^ | > 1 ) =AA = _1 2 

t— ¥OD t 2 
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Finally, our main theorem gives a strong law of large numbers for the catalytic BBM: 
Theorem 1.7 (SLLN). Let f : K — > R be some Borel-measurable bounded function. Then 

lim e-£* Y fiX?) = M x f f(x)pe-^dx P-a.s., 
ueN t JK 

where is the almost sure limit of the P -uniformly integrable additive martingale 

M t = J2 cxp{-/3|xn-i/? 2 i}- 

ueN t 

(Note: the martingale (M t )t>o will be discussed in detail in section^) 

One can observe that taking /(•) = 1 in Theorem 11.71 would give Lemma 11.11 and an even 
stronger result than given in Theorem ll.31 However, our proof of Theorem 1 1 . 71 relies on Theorem 
11.31 which in turn relies on Lemma 11.11 Thus, it will be necessary that we prove results in the 
order presented above. 

The rest of this article is arranged as follows. In Section 2 we recall some basic facts regarding 
the local times. We also introduce a Radon-Nikodym derivative that puts a drift towards the 
origin onto a Brownian motion. This will be useful in the subsequent analysis of the model. 
In Section 3, we recall some standard techniques for branching processes including spines and 
additive martingales. Section 4 is devoted to the proofs of Lemmas 11.11 and 11.21 We will prove 
Theorem 11.31 in Section 5, making use of the additive martingale (M t )t>o mentioned above. 
Section 6 contains the proofs of Theorem ll.41 Corollary II .51 and Lemma Tl. 61 Finally, in Section 
7 we give the proof of Theorem 11.71 this being largely based on extending the results found in 
Englander, Harris & Kyprianou |S]. 

2 Single-particle results 

Basic information about local times and the excursion theory can be found in many textbooks 
on Brownian motion (for example, see |14]V Also a good introduction is given in the paper of 
C. Rogers [TS]- Let us recall a few basic facts. 

Suppose (X t )t>o is a standard Brownian motion on some probability space under probability 
measure P. Let (L t )t>o be its local time at 0. Then (L t )t>o satisfies 

L ' = ^o2li 1 ^ 6 <- e ' e » ds 
for every t > 0. The next famous result is Tanaka's formula: 

\X t \= f sgn(X s )dX s + L t , 
Jo 

where 

1 if x > 
-1 ifa;<0 

In a non-rigorous way this can be thought of as Ito's formula applied to f(x) — \x\, where 
f'(x) — sgn(x), f"(x) = 2Sq(x) (where <5o is the Dirac delta function). Then one can think of L t 
as f Q 5 (X s )ds. 

Another useful result is the following theorem due to Levy. 



sgn{x) 
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Theorem 2.1 (Levy). Let (St)t>o be the running supremum of X . That is, St — sup 0<s<t X s . 
Then 

(S t , S t — X t )t>o — {L t , \X t \)t>o 



From Theroem 12.11 and the Reflection Principle it follows that Vt > 

L t = St = \X t \ = \N(0,t)\. 

It also follows that (Z t )t>o '■= {\Xt\ — L t )t>o — (Jq sgn(X s )dX s ) t >o is a standard Brownian 
motion under P, hence for any 7 G R 

exp {~/(\X t \ - L t ) - i 7 ^} = cxp {jZ t - i 7 2 *} , t > 
is a martingale. And more generally, for j(-) a smooth path we have the Girsanov martingale 



, r* 

Wt = exp • 



Tanaka ^ j ^ _ l _ j % ^ {s)ds } (2 .1) 

Used as the Radon-Nikodym derivative it puts the instantaneous drift sgn(X t )"/(t) on the process 
(X t )t>o- Let us restrict ourselves to the case 7 (-) = —7 < so that W puts the constant drift 7 
towards the origin on (X t )t>o- The following result is taken from [2]. 

Proposition 2.2. Let Q be the probability measure defined as 

^=^{-7(1^1-^)-^} , t>0, 

where (J r t)t>o *s the natural filtration of (X t )t>o- Then under Q, (X t )t>o has the transition 
density 



pit; a;, «) = — ■ exp 71 \x\ + \y\) 1 H — Lrtc[ = 

2y/2^i \ 2 It J 4 J \ V2l 

with respect to the speed measure 

m(dy) = 2e- 2 ^dy, 
so that 

Q x (X t e A) = [ p(t;x,y)m(dy). (2.2) 

J A 

Here Erfc(x) — -3= f°° e _u du ~ -^^e _;r as a; — > 00. 
It also has the stationary probability measure 

n{dx) = je^^dx. (2.3) 

3 Spines and additive martingales 
3.1 Spine setup 

In this section we give a brief overview of some main spine tools. The major reference for this 
section is the work of Hardy and Harris [10] where all the proofs and further references can be 
found. 
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We let (T t )t>o denote the natural filtration of our branching process as described in the 
introduction. We define T^ := a(Ut>oTt) as usual. 

Let us now extend our model by identifying an infinite line of descent which we refer to as 
the spine and which is chosen uniformly from all the possible lines of descent. It is defined in 
the following way. The initial particle of the branching process begins the spine. When it splits 
into two new particles, one of them is chosen with probability i to continue the spine. This goes 
on in the obvious way: whenever the particle currently in the spine splits, one of its children is 
chosen uniformly at random to continue the spine. 

The spine is denoted by £ = {0, £1, £2, • ■ ■ } ; where is the initial particle (both in the spine 
and in the entire branching process) and is the particle in the (n+l) st generation of the spine. 
Furthermore, at time t > we define: 

• node t {£,) := u e N t fl £ (such u is necessarily unique). That is, node t (t;) is the particle in 
the spine alive at time t. 

• m := \node t (S,)\. Thus n t is the number of fissions that have occured along the spine by 
time t. 

• £ t := Xf for u € N t n £. So (£t)t>o is the path of the spine. 

The next important step is to define a number of nitrations of our sample space, which contain 
different information about the process. 

Definition 3.1 (Filtrations). 

• Tt was defined earlier. It is the filtration which knows everything about the particles ' motion 
and their genealogy, but it knows nothing about the spine. 

• We also define Tt := o{Tt, nodet(£,)) ■ Thus T has all the information about the branching 
process and all the information about the spine. This will be the largest filtration. 

• Qt '■— °~{£,s ■ < s < t). This filtration only has information about the path of the spine 
process, but it can't tell which particle u € N t is the spine particle at time t. 

• Qt '■= o-(Qt, (node s (^) : < s < t)). This filtration knows everything about the spine 
including which particles make up the spine, but it doesn't know what is happening off the 
spine. 

Note that Gt C Qt C Tt and T t C T t . We shall use these filtrations to take various conditional 
expectations. 

We let P be the probability measure under which the branching process is defined together 
with the spine. Hence P = P]^. We shall write E for the expectation with respect to P. 

Under P the entire branching process (with the spine) can be described in the following way. 

• the initial particle (the spine) moves like a Brownian motion. 

• At instantaneous rate /3<5o(0 ^ splits into two new particles. 

• One of these particles (chosen uniformly at random) continues the spine. That is, it con- 
tinues moving as a Brownian motion and branching at rate (35o(-)- 

• The other particle initiates a new independent P-branching processes from the position of 
the split. 
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It is not hard to see that under P the spine's path (£t)t>a is itself a Brownian motion. We denote 
by (L t )t>o its local time at 0. 

Also, conditional on the path of the spine, (n t )t>o is a time-inhomogeneous Poisson process 
(or a Cox process) with instantaneous jump rate /3<$o(£t)- That is, conditional on Q t , k splits 
take place along the spine by time t with probability 

P(n t = k\g t )= { -^e-^. 

The next result (for example, see [10]) is very useful in computing expectations of various quan- 
tities 

Theorem 3.2 (Many-to-One Theorem). Let f(t) 6 mQ t . In other words, f(t) is Q t -measurable. 
Suppose it has the representation 



f(t) - ^ /«(*) 1 {™orfe t (e)=n}; 



ueN t 



where f u {t) G rnFt, the 



3.2 Martingales 



E { E fu{t))=E(f{t)e^ 

u£N t 



Since (£t)t>o is a standard Brownian motion we can define the following P- martingale with 
respect to the filtration (Gt)t>o using Proposition! 



Mf := e -m\+fiZ*-\lft , t >0. (3.1 
We also define the corresponding probability measure Qp as 

dQ/3 



dP 



= Mf , t > 0. (3.2) 

Qt 



Then under Qp, (£t)oo has drift (3 towards the origin and from Proposition 12.21 we know its 
exact transition density as well as its stationary distribution. 
Let us also define the martingale 

M t := 2 nt e~ f3Lt Mf ,t>0, 

which is the product of two P-martingales, the first of which doubles the branching rate along the 
spine, and the second puts the drift of magnitude /? towards the origin. If we define a probability 
measure Q as 

^| =Mt, t>0 (3.3) 

then under Q the branching process has the following description: 

• The initial particle (the spine) moves like a Brownian motion with drift /3 towards the 
origin. 

• When it is at position x it splits into two new particles at instantaneous rate 2/3<5o(x)- 



G 



• One of these particles (chosen uniformly at random) continues the spine. I.e. it continues 
moving as a biased random walk and branching at rate 2/35o (x). 

• The other particle initiates an unbiased branching process (as under P) from the position 
of the split. 

Note that although (|3.3p only defines Q on events in Ut>o^t, Caratheodory's extension theorem 
tells that Q has a unique extension on J 7 ^ := cr(Ut>o^Ft) and thus (|3.3p implicitly defines Q on 
J-oo. We then define Q :— so that 



dQ 
dP 



= M t 

Ft u£N t 



exp{(-/W|+L?->i)-/3L?} 
]T exp{-(3\X-\~\l3 2 t} , t>0. (3.4) 



uGN t 

{M t )t>o will be referred to as the additive martingale. 
3.3 Convergence properties of (M t ) t > 

The following theorem is a standard result for additive martingales in the study of branching 
processes. 

Theorem 3.3. (-M t ) t >o is P '-uniformly integrable and > P-almost surely. 

Proof. Recall the following measure-theoretic result, which gives Lebesgue's decomposition of Q 
into absolutely-continuous and singular parts. It can for example be found in the book of R. 
Durrett [B] (Section 4.3). 

Lemma 3.4. For events A e J 7 ^ 

Q(A)= / lim sup M t dP + Q{ A n { lim sup M t = oo } I . 

J A t—too ^ t— >oo ' 

Also a standard zero-one law, which can be found, for example, in |12j (see Lemma 3 and the 
proof of Theorem 2 that follows it) tells that P{M 00 > 0) G {0, 1}. Thus to prove Theorem 1331 
it is sufficient to show that 

lim sup M t < oo Q-a.s. (3-5) 

t— >oo 

Let us consider the spine decomposition of M t , another useful technique which can be found in 
EpfalGoo) =exp{-/3|&|-^ 2 <}+ J2 exp{-/?|e s J-^ 2 ^}, 

u<nodet (£) 

where we refer to the first term as spine(t) and the second term as sum(t). 

Recall that under Q, (£t)t>o is a Brownian Motion with drift (3 towards the origin and 
(|£t| — L t )t>o is a Brownian motion with drift —(3. Thus t -1 £t — > and t~ 1 L t — > (3 Q-a.s. Also 
spine(t) < 1 and 

sum(t)< e "^ S " <f>~^ 2s '\ (3-6) 

u<node t ((.) n=1 
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where S n is the n th birth on the spine. The birth process along the spine (rit)t>o conditional 
on the path of the spine is a time-inhomogeneous Poisson process (or a Cox process) with 
cummulative jump rate 2/3L t . Hence, Q-almost surely, n t ~ 2(3L t ~ 2(3 2 t, and so S n ~ (2/3 2 ) _1 n. 

Thus there exists some Q-a.s. finite random variable C > such that S n > Cn for all n. 
Substituting this into p.6[) we get 



oo 

1 t&t 



sum(t) < > e 3/3 Cn . 



»(*) < E 

n=l 

Therefore sum(t) is bounded by some Q-a.s. finite random variable. We deduce that 
lim sup E® I M t | Goo ) = lim sup ( spine(t) + sumit) j < oo Q-a.s. 

t— >oo V ' t— too V ' 

So by Fatou's lemma, Q-almost surely, 

( lim inf M t I ^ ) < lim inf f M t I ^ ) < lim sup E® ( M t I Go 



< OO. 



Then liminf t _ >00 M t < oo Q-a.s. and hence also Q-a.s. Since 1/M t is a positive Q-supermartingale, 
it must converge Q-a.s., hence 

lim sup Mt = lim inf Mt < oo Q-a.s. 

t-yoo t-taa 

completing the proof of the theorem. □ 

The next theorem is essential in the proof of the Strong Law of Large Numbers in the last 
section. 

Theorem 3.5. For p £ (1,2), (M t )t>o is L p - convergent. 

Proof. We use similar proof as found in [TU]. It is sufficient to show that E(M[) is bounded in t. 

e(m(\ = e(mp- x m^\ = e^^m?- 1 ) = eQ(m?-^ 
= EQ(EQ(M?- 1 \gJi)<EQ((EQ(M t \g ol ^ 1 ' 



by Jensen's inequality. Since for a, b > and q E (0, 1), (a + b) q < a q + b q , we see that 
r E^(Mt\g oo )Y~ = (spine(t) + sumit)) 1 '' 1 



< e -4r(p-i)t-<3(P-i)l« t | + E e " 



-(p-l)S„-/3(p-l)|«s„| 
uKnodet (£) 

And hence 



E(M^ < E®(e~'£ {p - 1)t - l3ip - 1)l ^ + E®( e-^ {p ~ 1)s ^ l3ip ~ 1) ^ s ^ 



u<nodet (£) 
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The first expectation is bounded by 1. The second one satisfies for e > small enough 
E °( J2 e -£(P-i)s.-0CP-i)l6 S «l) = £&( f 2j3L s e-^- 1 >-^- 1 ^ds) 

uKnoaet (£) 

^(2/3Z se -^(f- 1 )^(P- 1 'l^l)d S 



- ps ds 



ii 



< y £:Q^ e /3^ s -^-(p-i)s-/3(p-i)i« 3 |^ ds 

< I E Q^l e -PL s +l3\^\ + ^- Se (l+e)p(L s -\i s \)-^-ps^ ds 

3 2 (i+<e) 2 ,3 2 

e^-^ s "^ ps ds, 

which is bounded for e chosen sufficiently small (e < — 1). □ 

4 Expected population growth 

4.1 Asymptotic expected growth of |iV t | 
We prove Lemma 1 1 . 1 1 using the Many-to-One Theorem. 
Proof of Lemma From Theorem 13.21 we have 

25(M)=25(£l)=£(e^). 

ueN t 

Using the fact that L t = \N(0,t)\ it is then easy to check that 

where = P(iV(0, 1) < x). 

We can also find a good estimate of E{e^ Lt ) using the change of measure from (|3.2p . which 
is instructive for our purposes: 



Then, using the stationary measure, from (|2.3[) we have 
- / \ r°° r°° 



/3 2 t 



Thus 

□ 
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4.2 Asymptotic expected behaviour of 

Let us now prove that limt^oo \ \ogE{\N^\) = A A , where N t xt = {u E N t : X? > Xt} and 



A A = 



i/3 2 - /3A if A < 
-±A 2 if A > 

Proof of Lemma [T7B . Following the same steps as in the proof of Lemma \1 . 1 1 above we get 



4«t>At} 



ueN t 

= / e l3x p(t;0,x)m(dx)e^ 2t 
J\t 

At V2v^rt V t 2t ; 4 v ') 



where m(dx) and p were defined in Proposition [ 

Then for some functions ei(t) satisfying logej(i) = o(t) we have the following: 

f°° 1 -«i , , s -^t 
/ ■ e 2td£ = ei(t)e 2 , 

V2 7 Af v V2l y / \ £3We-^ + ^' +e 4 (t)e^* if A < /? 

where we have used that Erfc(x) ~ -^j^eT x as x — > oo and Erfc(ir) -)2asi^ — oo. Thus 



EQN?\) = 



e 6 (t)e-T* if A > /? 
eT^e-^^* if A < /3 



which proves the result after taking the logarithm and dividing by t. □ 

Again, we could have evaluated £^e /3 ^ t l{£ t>At }^ explicitly by using the joint density of £ t 
and L t (for example, see [12]): 

F(6 G ds, L ( G dy) = ^Lg exp { - ^ + y) * }dxdy , * G R, y > 0. (4.1) 

Lemmas II. II and 11.21 can also be proved via excursion theory (for example, see [2]). The proofs 
that we presented here (using the change of measure) in particular suggest the importance of 
the additive martingale (Mt)t>o in the study of the model. In the next section we shall see one 
simple application of this martingale. 
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5 Almost sure asymptotic growth of \N t \ 

In this section we prove Theorem 11.31 which says that log \N t \ ~ \0 2 t P-almost surely. 
Proof of Theorem \1.SX Let us first obtain the lower bound: 

liminf i^>l^ P _ a . s . (5.1) 

We observe that 

M t = ex P { ~ Wl - < l^tle"^ 2 *, 

hence logM t < log|iV t | - \fi 2 t and so ^logliN^ > ±/3 2 + t _1 logM t . Using the fact that 
lim.t_ J , 00 M t > P-a.s. from Theorem l3.3l we find that 

liminf^>V. 

t— s-oo t 2 



Let us now establish the upper bound: 



l0g|iV t | ^ 1 2 



limsup fal ' < -P z P-a.s. (5.2) 

t— >oo t 2 



We first prove (|5.2I) on integer (or other lattice) times. Take e > 0. Then 

P(\N t \e-W+* > c) < „ | e -« 

using the Markov inequality and Theorem ll.il So 

oo 

^P(|7V„| e -(^ 2 +^> e ) <0O . 
n=l 

Thus by the Borel-Cantelli lemma 

\N n \e- [ ^ f)2+e)n -> P-a.s. as n -> oo. 



Taking the logarithm we get 



Hence 



(-i/3 2 -e)n + log|Ay 



loglAU 1 9 

limsup 61 1 < -P 2 + e 



and taking the limit e — > we get the desired result. To get the convergence over any real- valued 
sequence we note that \N t \ is an increasing process and so 

logM < \t\ jog \ N [t] I 



* Til ' 
Hence 

logl^l , v log|iV rt1 | ^ 1 a 
limsup 1 < limsup < —p . 

t—yoo t t— too \t\ 2 

Combining (|5.2j) and (|5.1D now proves Theorem II .31 □ 



11 



6 Almost sure asymptotic behaviour of 1 7V^ A * | 



In this section we prove Theorem II .41 Namely, that 

t 2 

and 

\N^\ -> P-a.s. if A > ~. 

We break the proof into two parts. In subsection 6.1 we prove the upper bound and in subsection 
6.2 the lower bound. Also in subsections 6.3 and 6.4 present the proofs of Lemma [1.61 saying 
that limt^oo t~ 1 P(\N t xt \ > 1) = A A if A > §, and Corollary [T31 saying that lim^oo t^Rt = f , 
where R t is the position of the rightmost particle at time t. 



6.1 Upper bound 
Lemma 6.1. 



log|iV At | 
limsup 61 * 1 < A A P-a.s. 

t— >oo t 



The upper bound can be proved in a similar way to the upper bound on |iV t | (recall [5?2l . The 
main difference comes from the fact that (|7V t A '|) t >o is not an increasing process and so getting 
convergence along any real time sequence requires some extra work. 

Proof. Take e > and consider events 

An = { £ l { sup se[ „,„ + 1] X,>A, i} >e^ +£ )"}. 
u£N n+1 

If we can show that P(A n ) decays to exponentially fast then by the Borel-Cantelli Lemma we 
would have P(A n i.o.) = and that would be sufficient to get the result. 

By the Markov inequality and the Many-to-one theorem (Theorem 13. 2[) we have 



P(A n )<E( £ l {suPae[ „,„ +1] x ? > 



An} 



u£N r 



g -(A A +e)n 



+1 



= ^( e ^ + 1 l {sU p se[ „„ l+11 C,>A n} ) e " (A ^ e) " 
= ^(^" +1 l {5 „ +1+ C„>A n} ) e - (A ^ £) ", 

where :— sup sg [ n (£ s — is a sequence of i.i.d. random variables equal in distribution 

to sup sg [ 01 ] £ s and (£t)*>o is a standard Brownian motion under P. 

To give an upper bound on the expectation we split it according to whether |£ ra +i| is greater 
or less than (A — 5)(n + 1) for some small 5 > to be chosen later. 

^(^" + 1 %„ +1+ 6,>An})^ ( ^ +e) " 

+E[e^l {£n+1+in > a„}1{| ? „ +1 | < { x-8 ){n+ i )} )e- (A ^ n . (6.1) 
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Then from Theorem 11.21 we have 

Ilog(£(e^l u „ +1+fn > An} l m „ +1 | > ( A-,)(„ + i )} ) e - (Ai+e) " 
<ilog(^( e ^^l {|c „ +l| > (A _ 5)(n+1)} )e"( A ^^) 

= £ l0 S (2^(^ Z " +1 l {5n+1 > (A-^n+i)})) ~ (Aa + e) 
->A A _,5 - (A A + e). 

Since A A is continuous in A, A A _^ — (A A + e) < for 6 chosen small enough and hence the first 
expectation in (|6.ip decays exponentially fast. For the second expectation we have the following: 

^( e " £ " +11 {^ + 1+ ^ > A»}1{|^ +1 | < (A- 5 )(„ + l)}) e - (A ^ £)n 

> An}1{|?n+i| < (A _, )(n+1)} )e-(^^ 

<CE^ P (l {e „ +1 +f„ > Xn} 1 {|«„+i| < (A-5)(n+l)})e^" 

<^( 1 fe>^-A)})e K " 



=CQ^(6 > Sn+(6 -A)) 



where C = e^ 2 +( A " 5 ) and if = \f3 2 + f3{\ - 6) - (A A + e). However Q^i > <5n + (S - A)) 
decays faster than exponentially in n. To see this observe that for any 9 arbitrarily large 

where 

< (e^f(e^f <oo 

using the Cauchy-Schwarz inequality and the fact that L\ = £i = \N(0, 1)| under P. Thus we 
have shown that the expectation in (16.11) and consequently P(A n ) decay exponentially fast. 
So by the Borel-Cantelli lemma P(A n i.o.) = and P{A„ ev.) = 1. That is, 

£ 1 {su Pse[n , n+1] x- > An} < e^ +e >™ eventually. 

u£N n+1 

So there exists a P-almost surely finite time T e such that Vn > T e 

E l{su Pse[ „,„ +1] x ? >A n} <e( A ^)". 



Then 



|^V t At | < 2J 1 {sup 3e[L tj, uj+h*? ^ A L*J> 
.j^Aij < e (A x +e)L*J fort>T e + l, 
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which proves that 



log|7V t At| 



limsup < Aa P-a.s. 



□ 

Remark 6.2. Since \N xt \ takes only integer values we see that for X > ^ the inequality 

limsup Sl 1 1 < A A < 

t— >oo I 

actually implies that \N xt \ — > P-a.s. 



6.2 Lower bound 

Before we present the proof of the lower bound of Theorem [Lj] let us give a heuristic argument, 
which this proof will be based upon. 

Take A > 0. Suppose we are given some large time t and we want to estimate the number of 
particles u G N t such that > At. 

Let p £ [0, 1]. At time pt the number of particles in the system is \N pt \ ss e^ 13 pt by Theorem 
11.31 If we ignore any branching that takes place in the time interval (pt, t] then each of these 

particles will end up in the region (-co, —Xt] U [Xt, oo) at time t with probability > e 2 < 1 -p) 
using the standard estimate of the tail distribution of a normal random variable. 

Thus a crude estimate gives us that the number of particles at time t in the region (— oo, — Xt]U 
[Xt, oo) is 

> e -3rT=57* x |_zv pt | « e-air^T' x e^ 2pt . 
The value of p which maximises this expression is 



P 

and then 



log(e w-?)' x \Npt\j 



if A > p 
j ifA</3 



-±A 2 if A > /? _ 
p=p . ~ \ \P 2 -pX if A < p 



t 

Let us now use this idea to give a formal proof of the following lemma. 
Lemma 6.3. TakeX<§. Then 

t-^oo t 2 

Proof. Take p := 1 — G (|, 1). For integer times n we shall consider particles alive at time pn 
(that is, particles in the set N pn ). 

For each particle u £ N pn we can choose one descendant alive at time n + 1. Let iV n +i be a 
set of such descendants (so that |-/V„ + i| = |iV pn |). 

Then, for u G N n +i, paths (Xl l ) te ^ pn correspond to independent Brownian motions (started 
at some unknown positions at time pn). Note that, wherever particle u is at time pn, 

P(\X?\ > Xs Vs G [n, n + l})> er^) n = e ~^ Xn =: q n (X) 
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using the tail estimate of the normal distribution. Take a small 6 > to be specified later. Then 
by Theorem II .31 

\N n+1 \ = |iVp„| > e ( ^ 2 P~ 5 ^ eventually. 
To prove Lemma 16.31 we fix an arbitrary e > and consider the events 

Bn := { ^2 1 {\ x l , \> Xs Vse[n,n+i]} < e (AA ~ e) ™j. 

u£N n+1 

We wish to show that P(B n i.o.) — 0. Now, 
p(B n n{\N n+1 \ >e ^- s >]) 
=p({ Y, M\x?\>Xs v se[ n,„+i]} < e (AA " e)n } n {\N n+1 \ > e^- 6 ^}) 



<p( U<e (A ^ e)n 
where A^s are independent events with P(Ai) > q n (X) Vi . Then 



P J2 l Ai <e^-^ n )=P(e-^ 



i=i 



> e 



= e e<A ^ e> " J] 25 ( 
i=i 

< e e^-^ n ( 1 _p (Ai)(l _ e -l ) 

<e e(i """n(l-9»W(l-^ 

= cxp j e ( A ^)« _ (l _ e - 1 )g„(A)e^> 92 P- 5 )«} 
= exp | e ( A A-^)" _ (l _ e -i) e (A A -«)»\ 
This expression decays fast enough if we take 5 < e. Thus 

p(p n n {|iv n+1 | > e (^ 2 p- 5 )™} i.o.) = o. 

And since p({|iV n+ i| > e ih0 2 P- s ) n ) ev. ) = 1, we get that p(p„ i.o. ) = 0. That is, 
Y l{|X«|>As Vse[n,n+i]} > e (AA ~ e)n for n large enough P-almost surely. 



+ 1 



Now, since the process is symmetric, the probability that a particle u E N n +i such that > 
As Vs € [ri, n + 1] actually satisfies X™ > As Vs <E [n, n + 1] is |. So applying the usual Borel- 
Cantelli argument once again we can for example prove that for some constant C € (0, i) 

l{x a ">As Vse[n,n+i]} > Ce^ -6 '" for n large enough P-almost surely. 



U£N T 



+1 
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Then for t large enough 



|JV*| = J2 1 W>At} > E Va 6[L *j,L*j+i]} > Ce^-)W. 



Thus 

liminf«n>A A . 



Lemmas 16.11 and 16.31 together prove Theorem 11.41 
6.3 Decay of P(\N^\ > 1) in the case A > f 

Theorem [L4] told us that if A > § then |iV A *| — > 0. Let us also prove that in this case 



□ 



log P(\N t xt \> 1) = _1 a2 
i A 2 

Proof of Lemma[TE Trivially P(|iV A *| > 1) < £|iV t At |. Hence by TheoremO 

logP(|iV t A '| > 1) 

t— >oo £ 

For the lower bound we use the same idea as in Lemma 16.31 Let us take 



if A > j3 
l ' ~ ^ 1 - | if A < /3 

We define a set N t as in Subsection 6.2 That is, for each particle u € N pt we choose one descendent 
alive at time t (so that Nt C iVj, |iVf| = |AT pt |). Then for each u £ N t wherever it is at time pt 
we have 

PQX?\ > Xt) > e"^* =: Pt (X). 

Then 

p(|iV At |>l)>ip(|7V t ±A *|>l), 
where N^ xt := {u e N t : \X?\ > Xt}. Thus for a small 5 > to be specified later we have 



P(\N t xt \ > l) > ip(|iV t ±At | > 1, \N pt \>eWP- s ») 
>\p{ (J i\ x t\>^}, \N pt \>n t (6) 

u£N t 

1 -(l-(l- Pt (X)r {S) )p(\N pt \>n t (5)) 



~ 2 

By Theorem [L3l P(\N pt \ > n t (S)) -> 1, so we can ignore this term. Then 

(l - (1 - Pt (A))^) = n t (S) Pt (X) ~ ( n f )p t (A) 2 + ( n f)^(A) 3 -- 

> n t (S) Pt (X) - n t (S) 2 Pt (X) 2 (l + n t (S) Pt (X) + n t (S) 2 Pt (X) 2 + 
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Note that for 8 small enough 



n t (8)p t (\) = eWr-We-itt' = e (A ^ 5)t « 1. 
Hence P(\N t xt \ > 1) > (^P(\N pt \ > n^S^e^-^ 1 + o(e^- 5 ^) and therefore 

liminf logP(l^l<D^ AA 

t— >oo t 

This completes the proof of Lemma 11.61 □ 
6.4 The rightmost particle 

Observe that the number of particles above the line At grows exponentially if A < 4 and is 

— I [run j » -i*. \i i 'i ill ill \ ill I Iw.iw/.vn I I w<< 07. i ili'ii 



eventually if A > ^. Hence, as a corollary of Theorem II .41 we get that 



> — P-a.s., 

t 2 

where (Rt)t>o is the rightmost particle of the branching process. 

Proof of Corollary IT731 Take any A < f. By Theorem H~4l |7V t At | > 1 Vi large enough, so R t > Xt 
for t large enough. Thus liminft^oo i _1 P t > A P-a.s. Letting A ^ we get 

hminf — > — P-a.s. 

Similarly, if we take A > § then by Theorem II .41 |jV f At | = Vt large enough and so i?t < A£ for t 
large enough. Hence limsup^^ i _1 i£t < A P-a.s. So, letting A \ ^ we get 

r R * <- & p 
hmsup — < — P-a.s. 

t— foo c 2 

and this proves Corollary 1 1.5 1 □ 

Note that the rightmost particle (that is, the extremal particle) in our model behaves very 
differently from the rightmost particle in the model with homogeneous branching. 

In the BBM model with homogeneous branching rate (3 there is a particle staying near the 
critical line \f2fit all the time. (The word particle here is a bit ambiguous since we are really 
talking about an infinite line of descent, but this is a common description.) 

On the other hand in the BBM model with branching rate f3So(x), since branching only takes 
place at the origin, no particle can stay close to any straight line At, A > all the time. The 
optimal way for some particle to reach the critical line at time T is to wait near the origin 
until the time ^ in order to give birth to as many particles as possible, and then at time % one 

of ss eT T particles will have a good chance of reaching at time T. 
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7 Strong law of large numbers 

Recall the additive martingale M t = e~~ l Y^ueN e ~' 9 ' X "'j t > from (13.41) and the measure 
7r(dx) = f3e~ 2/3 \ x \dx from Proposition 12.21 In this section, we shall discuss Theorem 11.71 which 
says that for a measurable bounded function /(•) 

a 2 f°° 

J -™ 

/CO 
f[x)e^Tr(dx) P-a.s. (7.1) 
-oo 

Observe that convergence in L 1 is trivial by the Many-to-One theorem and the uniform integra- 
bility of (M t ) t > : 



u£N t 

= P>(/(6)^ l?t| ) f f(x)e^n(dx). 



-P\£t\+PL t -2- 



Also the Weak Law of Large Numbers for this model have been proved by J. Englander and D. 
Turaev in [7] . In particular they have given the law of . 

As a simple corollary of Theorem 1 1.71 we get by taking /(•) = 1 that 



N t \e-^ 2t -» 2Afoo P-a.s. (7.2) 



This should be compared with results in Lemma 11.11 and Theorem 11.31 
Dividing (|7.1[) by (17.21) we get an alternative form of Theorem 11.71 

S " e ^ m "> ~ / me^Mdx) = U f(x)e-^dx P-a.s. 



The Strong Law of Large Numbers was proved in [5] for a large class of general diffusion processes 
and branching rates fi(x). In our case the branching rate is a generalised function /38o(x), which 
doesn't satisfy the conditions of [8] . Nevertheless we can adapt the proof to our model if we take 
the generalised principal eigenvalue A c — 4j- and the eigenfunctions (j>(x) — e~^ x \ 4>(x) = /Se - ^' 
in [8]. Also the proof relies on the LP convergence of the martingale (M t ) t >o and the linear 
asymptotic growth of the rightmost particle which we have derived earlier in this article. 
Let us now finish the article with the proof of Theorem 11.71 

Proof of Theorem \1.7\ Take B a measurable set. As it will be shown later, it is sufficient to 
prove the theorem for functions of the form f(x) = e~^' x 'l^ xe B}- For such set B let 

u£N t ueN t 

So if B = R then we would have Ut = M t and generally Ut < M t . We wish to show that 
U t tt(P)M 00 ( = / f(x)e^Tr(dx)M 00 ^ as t oo. 
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The proof can be split into three parts. 
Part I : 

Let us take K > 0. At this stage it doesn't matter what K is, but in Part II of the proof we 
shall choose an appropriate value for it. Let m n := Kn (using the same notation as in [8]). Also 
fix 5 > 0. We first want to prove that 



lim 



U, 



)6-e(u { 



= P-a.s. 



We begin with the observation that 



vm > o u s+t = 

ueN t 



tjj(u) 
s ' 



(7.3) 



(7.4) 



where conditional on Ft, Us are independent copies of U s started from the positions Xf. 
To prove (|7.3p using the Borel-Cantelli lemma we need to show that for all e > 

oo 

^2 P {\ U (n+m n )8 - E(u {n+mn) s\F n s^ > ej < co. (7.5) 

n=l 

Let us take any p £ (1, 2). Then 

Next we shall apply the following inequality, which was used in the proof of the SLLN in [5] and 
can also be found in [I]: if p £ (1, 2) and Xi are independent random variables with KXi = (or 
they are martingale differences), then 



(7.6) 



i=l 



Then by E3 



Us+t ~ E(u s+t \F t ) = J2 e-^(t/ s (u) - E(U^\F t 



ueN t 



where conditional on Ft, U^ — E(Us U ^\Ft) are independent with mean. Thus applying (|7.6p 
and Jensen's inequality we get 



E(\u s+t -E(u s+t \F t ) 



Vt) <2Pe-^ 4 E(\U^ -E(U^\F t )\ P \F t 

ueN t 

<2P e -p£* E(2 p - 1 (\U^\ p +\E(UM\F t )\ P )\F t S ) 



<& e -v4t J2 E^^U^ + E(\U^\ p \F t ))\J 7 t 

ueN t 

= #P e -p£* £ E (\U^\ p \F t ). 



(7.7) 



ueN t 



19 



Hence by (|777|) 

OO 



M„ 



n=l u£N Sn 
= 2*£V*^ n £( £ e -*l- x ?»lB (Jlf mn4 ) P 

n=l ueN Sn 

00 2 2 



n=l 



where C is some positive constant and we have applied the Many-to-One Theorem (Theorem 
13. 2j) and and Theorem 13.51 in the last inequality. Since p > 1 the sum is < 00. This finishes the 
proof of (|7.5[) and hence (|7.3I) . 
Part II: 



Let us now prove that 



lim 



f^(n+m„)5l-^n«) - 7r(5)M 



= P-a.s. 



(7.8) 



Together with (|T.3[) this will complete the proof of Theorem 11.71 along lattice times for functions 
f(x) of the form e~ /3 ' a: 'l{ xe B}. 
We begin by noting that 



uGN t 



= ]T e~^E x "U s 

u€N t 

ugff, «GiV s 

uefft 

: ^ e 4M|xn £ gf( 1{fjefl} ) 

«GAf t 



e 2 



4-t-/3|X t "| 



p(s,X?,y)m(dy), 



where Qp and p(-) were defined in (|3.2p and Proposition 12.21 Thus 

S(c/ (n +m»)5|^n«) = ]T e'^ n5 '^ x ^ f p(m n S, X^y)m(dy). 



(7.9) 



uGiV„<5 

Recalling that m„ = ifn where K > we have 



^(c/(„ +m „) 5 |J-„ 5 ) = ^ e"^" 5 "' 3 !^! / P (KnS,X: s ,y)m(dy). 



ueN„ s 
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Now choose M > § and consider events 



C n := [\Xl 5 \ < Mn6 Vu e iY n(5 }. 



Then 



E e-^ n5 -^ |X « sl / p(Kn5,X^,y)m(dy) 



ueiV n6 " a 

+ e-^ n6 ~^\ [ p(KnS,X^, y )m(dy)l Cn . 

The first sum is for n large enough by Corollary 1 1.5 1 (or even earlier by Theorem II .4[) . To deal 
with the second sum we substitute the known transition density p(-): 



p(KnS,X% s ,y)m(dy) l Cn 



1 



K^J\ X ns\ + \v\-PKn6 



2 ErfC 



2 2An<5 
2/3,5,1 dy l c „. 



Then for any given M > § we can choose K > and hence 

< exp | (pM - yifj njj x C' -^Oasn-> oo, 
where C" is some positive constant and 



2Kn8 



Ay l c 



^ e -2^|w| dy = 7r (B) asrwoo 



since Erfc(x) -> 2 as i -> — oo and lc„ -> 1 as n -> oo. Then going back to (|7.9[) and we see 
that 



and so also 



lim 



lim 

n— >oo 



E[U( n+mn ) S \F n s) -n(B)M nS 



E 



= P-a.s. 



P-a.s. 



As it was mentioned earlier parts I and II together complete the proof of Theorem 11.71 along 
lattice times for functions of the form f(x) = e~^ x ^ls{x). To see this put together (17.3[) and 
(EU) to get that 



That is, 



lim 



lim 

n— f oo 



U(n+m n )8 - Tt(B)M 

U n (K+i)s - n(B)M a 



= P-a.s. 



P-a.s. 



21 



Then K + 1 can be absorbed into 6 which stayed arbitrary throughout the proof. Also as it was 
mentioned earlier we can easily replace functions of the form e~^lB(x) with any measurable 
functions. To see this we note that given any meausurable set A and t\ > we can find constants 
c 1 , . . . , c n , Ci, . . . , c n and measurable sets A\, . . . , A n such that 

n n 

(^QUWe-* 1 ) -ei < l A (x) < J2^AMe~ m 
i=l i=l 

and 

n n 
i=l i=l 

Similarly given any positive bounded measurable function / and e 2 > we can find simple 
functions / and / such that 

/» - e 2 < f(x) < f(x) 

and 

1(X) < f(x) < fjx) + 6 2 . 

Thus given any positive bounded measurable function / and e > we can find functions f e {x) 
and f t (x), which are linear combinations of functions of the form e _/3 l x ll^(a;) such that 



and 
Then 



f £ (x)<f(x)<r(x)+e. 



/oo poo 
f(x)f3e~ f3lxl dx < / (f(x) + e)(3e- f3 ^dx + 2e 
-co J —oo 



and hence P-almost surely we have 

limsupe"^" 5 Yl /(^)<linisup e -^ n5 £ fffi) 



-CO 



< 



/CO 
f^e^^lda; 
-CO 

Moo( / /(x)/3e-' 3 l a: ldx + 2e). 



Since e is arbitrary we get 

2 /*00 

limsupe--" 5 V /(X^KMn f(x)p e -^dx. 

n— s-oo J— oo 



MSiV„,s 

Similarly 



2 r°° 
liminfe"^" 5 V /«) > / /(x)/3e~^ldx. 
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Also any bounded measurable function / can be written as a difference of two positive bounded 
measurable functions. This completes the proof of Theorem 11.71 with the limit taken along lattice 
times. Now let us finish the proof of the theorem by extending it to the continuous-time limit. 
Part III : 

As in the previous parts of the proof it is sufficient to consider functions of the form f(x) = 
e~P\ x \l B (x) for measurable sets B. 

Let us now take e > and define the following set 



B'l.r) , ZM ( -|z|-ilog(l + e), | x | + Ilog(l + e ; 



Note that y E B e (x) iff y E B and > ^j-^r- Furthermore, for 5, e > let 



-b( x ) 



L {X,"6B f (i) VsG[0,<5] VuGiVa} 



and 



Then for t E [nS, (n + 1)5} 



U t = 



t = e 



: E< 

ueN t 



-p\x?\ 



= Y, e-^UpL > E U i-ns ^b(Ks) 

P 2 r » fl 2 r p-^l^nsl 



1 + e 



(7.10) 



because at time t there is at least one descendent of each particle alive at time nS. Let us consider 
the sum 



e 2 



n S 



E 



u£N„ 



Note that 



^^(Xnfi) are independent conditional on J- n s, 



e -?\K s \^{Xl 5 ) 



F n s = e 



j2 e -«idrra, 



and 



hm e~^ nS £ e^ x ^ S B{Ks) = [ &(xM^)M 

n— yoo * — ■* / 



(7.11) 
(7.12) 

(7.13) 



u£N nS 



The last equation follows from the SLLN along lattice times which we already proved. Also we 
should point out that if we further let S — > 0, i% e {x) will converge to 1 B (%) and (|7.13p will 
converge to tt(B)M 00 . Our next step then is to show that 



lim 

n— >oo 



E * 

u£N nS 



-i iS 



E 



u£N n 



(7.14) 
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In view of (|7.11[) and (|7.12l) we prove this using the method of Part I. That is, we exploit the 
Borel-Cantelli Lemma and in order to do that we need to show that for some pe (1,2) 



u£N„ 



B K^nS, 



°° / 2 ,2 

n=l ^ ueN„e uEN n 

A similar argument to the one used in Part I (see (|7.7p gives us that 

00 / 2 2 



< 00. 



71=1 
OO 



n=l 



where S^ e (X% s ) is an indicator function and therefore raising it to the power p leaves it un- 
changed. Using once again the Many-to-One Lemma and the usual change of measure we get 

OO 2 

Y&e-^E^ £ e-*l*ai#(X£)) 

n=l ueN nS 
00 2 



u£N„ 



71=1 

OO 



u£N n 



^2 2 2p e~ {p ~ 1) ^ nS E^f( K Y e"' 9(p " 1)|x "' sl 



< 00. 



71=1 ueN nS 
Thus we have proved (|7.14p . which together with f)T. 13[) implies that 

liminfe"^ £ B -^H*-(J%) = Kminf e"^ £ e^^l^TO 



^ £ (a;)7r(da;)M 00 



Putting this into (|7.10[) and letting n = ||J gives us 



e 2 

lim inf C/j > 

t=>oo 1 + e 



^ e (x)7r(dz)M 



Letting 5, e \ we get liminf t _ i . 00 C7< > 7r(_B)Af 00 . Since the same result also holds for B c we 
can easily see that limsup^,^ Ut < n(B)M 00 . Thus 

lim U t =ic(B)M 00 . 

Then the same argument as at the end of Part II of the proof extends the result for functions of 
the form lB(x)e~P' x ' to all bounded measurable functions. □ 
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